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Abstract:  The  collection  of  light  at  very  high  numerical  aperture  allows 
detection  of  evanescent  waves  above  the  critical  angle  of  total  internal 
reflection  in  solid  immersion  lens  microscopy.  We  investigate  the  effect  of 
such  evanescent  modes,  so-called  forbidden  light,  on  the  far-field  imaging 
properties  of  an  aplanatic  solid  immersion  microscope  by  developing  a 
dyadic  Green’s  function  formalism  in  the  context  of  subsurface 
semiconductor  integrated  circuit  imaging.  We  demonstrate  that  the 
collection  of  forbidden  light  allows  for  sub-diffraction  spatial  resolution  and 
substantial  enhancement  of  photon  collection  efficiency  albeit  inducing 
wave-front  discontinuities  and  aberrations. 
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1.  Introduction 

The  growing  demand  for  high  spatial  resolution  in  optical  inspection  for  failure  analysis  of 
semiconductor  integrated  circuits  (ICs)  has  aroused  interest  in  employing  aplanatic  solid 
immersions  lenses  (aSILs)  [1-7].  In  such  applications,  aSILs  are  placed  in  intimate 
mechanical  contact  with  the  polished  back-side  of  an  IC  chip,  hence  transforming  the  silicon 
substrate  into  a  high  refractive  index  immersion  medium  (nSi  =  ~3.5)  for  high  numerical 
aperture  (NA)  sub-surface  imaging. 

The  circuit  features  fabricated  on  the  front-side  of  the  silicon  substrate  such  as  resistors, 
capacitors,  interconnects  and  transistor  gates  are  surrounded  by  an  insulating  medium  with  a 
lower  refractive  index  (typically,  nins  =  ~1.5)  compared  to  the  silicon  immersion  index.  The 
electromagnetic  boundary  conditions  at  the  dielectric  interface  between  the  silicon  substrate 
and  insulating  media  impose  that  the  evanescent  waves  originated  from  the  objects  in  the 
insulating  medium  are  transformed  into  propagating  waves  in  silicon  immersion  medium  at 
angles  higher  than  the  critical  angle  of  total  internal  reflection  (TIR)  [8].  Such  propagating 
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waves,  so-called  forbidden  light,  can  be  collected  by  the  aSIL  and  thus  contribute  to  the  far- 
field  imaging  [9,  10].  In  modern  IC  manufacturing  technology,  most  important  circuit  features 
reside  within  the  depth  of  the  decay  length  of  evanescent  waves  therefore  the  forbidden  light 
has  important  ramifications  for  high  NA  imaging  of  ICs  [11,  12].  Recently,  the  imaging 
characteristics  of  aSIL  based  microscopes  have  been  theoretically  investigated  [13-16].  Due 
to  their  limited  scope  to  imaging  objects  buried  in  a  homogenous  immersion  medium,  these 
previous  analyses  are  not  directly  applicable  for  assessing  the  high-NA  imaging  performance 
of  aSIL  through  stratified  dielectric  media  such  as  semiconductor  IC  chips. 

In  this  article,  we  investigate  the  imaging  performance  of  a  high-NA  aSIL  microscope 
through  developing  a  dyadic  Green’s  function  formalism  for  the  analysis  of  the  forbidden 
light.  It  is  shown  that  the  forbidden  light  allows  sub-diffraction  spatial  resolution  [17]  and 
approximately  a  factor  of  5  improvement  in  the  photon  collection  efficiency  for  imaging  the 
metal  levels  closest  to  silicon  as  compared  to  those  at  a  distance  of  more  than  a  half¬ 
wavelength  into  the  insulating  medium.  Furthermore,  the  evanescent  nature  of  the  collected 
light  at  supercritical  angles  leads  to  wave-front  aberrations.  Investigation  of  imaging 
performance  in  the  vicinity  of  the  dielectric  interface  through  an  electromagnetic  model  is 
critical  to  assess  the  limitations  of  the  aSIL  microscopy  for  semiconductor  failure  analysis  and 
is  applicable  to  imaging  in  quantum  optics  [18],  biophotonics  [19]  and  metrology  [20]  as  well, 
as  pure  ray-tracing  models  cannot  account  for  the  behavior  of  the  evanescent  waves  discussed 
here  [21,  22], 

2.  Theory 

The  theoretical  analysis  is  divided  into  two  steps.  We  first  develop  a  model  of  the  angular 
spectrum  of  the  radiation  from  a  dipolar  object  located  in  the  vicinity  of  the  aplanatic  point  of 
the  aSIL.  Secondly,  we  present  the  details  of  the  dyadic  Green’s  function  of  an  aSIL 
microscope  to  investigate  the  high-NA  far- field  imaging  of  buried  objects  beyond  the 
interface  between  the  substrate  and  the  insulating  media. 

The  schematic  illustration  of  the  problem  is  shown  in  Fig.  1.  A  dipolar  object  with  an 

electrical  dipole  moment  ji  =  (jix,jly, jlz)  is  located  at  |rd |  =  {yxd,yd,zd)  in  the  insulating 

medium  that  has  a  lower  refractive  index  (nim)  than  the  silicon  immersion  medium  ( nasn )•  The 
propagating  waves  in  the  +  z  direction  of  the  coordinate  system  of  the  microscope  undergo 
refraction  at  the  planar  interface  between  the  immersion  and  insulating  media  and  are 
channeled  towards  the  subcritical  angle  cone  upon  refraction  (allowed  light,  green  in  Fig.  1). 
The  evanescent  waves  on  the  other  hand,  are  transformed  into  propagating  waves  beyond  the 
critical  angle  at  the  interface  (forbidden  light,  red  in  Fig.  1).  The  outward  propagating  fields  in 
the  immersion  lens  medium  comprised  of  allowed  and  forbidden  light  are  then  refracted  on 
the  spherical  surface  of  the  aSIL  and  collected  by  the  backing  objective  lens.  The 
conventional  ray  tracing  models  cannot  account  such  far-field  characteristics  originating  from 
non-propagating  modes  and  thus  an  electromagnetic  solution  is  strictly  required  for  the 
problem.  The  electromagnetic  solution  to  the  problem  illustrated  in  Fig.  1,  can  be  formulated 
using  the  Green’s  function  formalism.  Assume  that  the  dyadic  Green’s  function  G0(r,rd) 

defines  the  electric  field  E(r)  of  an  electric  dipole  jl  located  at  rd  in  the  homogenous 
insulating  medium: 

E(r )  =  CO2  /J0G  0  (r  ,rd )  •  jl  (1) 

where  co  and  /./„  represent  the  oscillation  frequency  of  the  light  and  vacuum  permeability, 
respectively.  The  Green’s  function  satisfying  the  homogenous  space  solution  is  given  as  [23]: 
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where  kins  refers  to  the  wave-number  ( In nim  /  A0 )  on  which  the  subscript  shows  the 
component  in  a  given  coordinate  axis.  Equation  (2)  allows  us  to  express  the  fields  of  a  dipole 
in  Cartesian  coordinates  in  terms  of  plane  and  evanescent  waves  in  z  >  zd  assuming  the 
insulating  medium  is  homogenous.  The  next  step  is  modifying  the  Green’s  function  to 
account  for  the  planar  dielectric  interface  and  retain  the  far-field  terms  to  simplify  the  solution 
in  Eq.  (2)  following  the  procedure  in  Ref  [23]: 
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The  Green’s  function  in  Eq.  (3)  allows  deriving  the  far-fields  in  spherical  coordinates  at  a 
given  location  in  the  immersion  medium  ( r  »  A„  and  r  »  |rd  | ).  The  terms  t'"s  and  t‘"s  are 

the  Fresnel  coefficients  for  transmission  at  the  planar  interface  between  the  immersion  and 
insulating  media.  The  remaining  task  in  the  problem  is  accounting  the  refraction  on  the 
spherical  surface  of  the  aSIL  in  order  to  find  the  electric  field  on  the  Gaussian  reference 
sphere.  Following  Ref  [16],  the  electric  field  on  the  Gaussian  reference  sphere  is  obtained: 
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where  &o/yand  kim  are  the  wave-numbers  of  the  light  in  the  objective  and  insulating  media, 
respectively;  fobJ  is  the  focal  length  of  the  objective,  kzim  and  kzaSIL  are  the  longitudinal 
components  of  the  wave-vectors  in  the  insulating  and  immersion  media,  respectively;  d  is  the 
distance  of  dielectric  interface  from  the  aplanatic  point;  9aSIL  and  6obj  represent  the  polar 
angles  with  respect  to  the  aSIL  and  the  objective  coordinate  centers,  respectively. 
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typical  aSIL  based  imaging  system,  the  light  collected  by  the  backing  objective  is 
on  a  light  detector  such  as  a  camera  by  a  tube  lens.  Assuming  the  aplanatic  conditions 
active  and  the  tube  lenses,  the  dyadic  Green’s  function  describing  the  fields  on  the 
due  to  an  arbitrarily  oriented  point  dipole  source  near  an  interface  is  obtained  below: 
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where 


P  =  \]x2+y2 ;  (p  =  tan-1  (y/x) 

*  =  -  (^de,  sin  #det*de,  +  sin  daSILXd  ) 
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The  derivation  steps  between  Eqs.  (4)  and  (5)  are  skipped  for  the  sake  of  brevity;  and  these 
details  can  be  found  in  Ref  [16].  In  Eqs.  (5)  and  (6),  kia  is  the  wave  number  in  the  detector 

space,  fdet  is  the  focal  length  of  the  tube  lens  (/*, ,  »f0bj),  >idet  is  the  refractive  index  of  the 
detector  medium,  rdet  =  (xdet,  ydet,  zdet)  is  the  location  on  the  detector  plane.  J(p)  refers  to  the 
Bessel  functions  of  a  given  order.  The  relationship  between  the  polar  angles,  9aSIL  ,  9obj  and 

9ia ,  is  as  follows: 


nobj  sin  9aSIL  —  naSIL  sin  9obj 


/det  sin  9det  =  fobj  sin  9C 


^ obj 


(7) 


>hn,  S'n  dlns  =  naS,L  sin  9aSI,_ 


The  Green’s  flinction  expressed  in  the  cylindrical  coordinates  of  the  detector  in  Eq.  (5). 
allows  us  to  obtain  the  far-field  image  of  an  arbitrarily  oriented  dipole  located  in  the 
proximity  of  the  aplanatic  point  of  the  aSIL  near  the  planar  interface.  Note  that  the 
formulation  of  the  Green’s  function  is  valid  provided  that  the  radius  of  the  aSIL  is  much 
greater  than  the  wavelength  of  light  and  the  distance  of  the  dipole  location  from  the  aplanatic 
location  ( RaSII  »  A,,  and  RaSIL  »|?d|)  [16],  typically  satisfied  in  most  practical  applications 


[5-7,  22,  24-26].  A  more  sophisticated  numerical  model  is  required  for  cases  in  which  such 
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assumption  is  not  valid  [27,  28].  The  objective  and  the  tube  lenses  are  assumed  to  be  ideal 
meaning  that  they  are  free  of  optical  aberrations  and  have  a  transmittance  of  unity  in 
accordance  with  the  previous  studies  [14-16].  Such  practical  details  can  be  incorporated  into 
the  model  provided  that  the  empirical  data  is  available.  Furthermore,  the  developed 
electromagnetic  model  is  intended  to  solve  the  light  propagation  from  object  space  to  the 
detector.  The  stray  light  originating  from  the  reflections  from  lens  surfaces  or  chip  package 
are  not  considered  in  the  image  formation  process  in  accordance  with  the  previous  studies  [3- 
7,  14-16], 


R  naS|L/nobj 


Fig.  1.  Schematic  of  the  problem.  GRS  and  R  stand  for  Gaussian  reference  sphere  of  the 
objective  and  radius  of  the  aSIL,  respectively.  Green  and  red  zones  denote  the  allowed  and 
forbidden  light  regions,  respectively. 


3.  Results  and  discussion 

In  the  following,  we  first  elucidate  the  characteristics  of  the  forbidden  light  through  studying 
the  field  amplitude  and  phase  profile  on  the  objective  pupil.  Secondly,  we  investigate  the 
point  spread  functions  to  understand  the  dependence  of  spatial  resolution  and  light  collection 
efficiency  on  the  depth  of  the  objects.  A  case  of  two  non-interacting-dipoles  is  briefly 
discussed  to  include  the  effect  of  coherence  in  the  spatial  resolution.  Finally,  we  discuss  the 
simulated  image  of  two-dimensional  objects  buried  in  the  insulating  medium  such  as  the 
metal  wiring  in  IC  chips. 

Following  parameters  are  used  in  all  simulations:  nasiL  =  3.5,  nins  =  1.53,  nabj  =  \,fobj  = 
10mm,  Rasn  =  2.35mm,  the  free  space  wavelength  ).n  =  1340nm  and  objective  numerical 
aperture  NA0bj  =  0.278.  These  values  reflect  typical  experimental  parameters.  The  objects  are 
situated  on  the  same  transverse  plane  of  aplanatic  point  of  the  aSIL  (zd  =  0),  regardless  of  the 
their  axial  distance  from  the  dielectric  interface  unless  it  is  stated  otherwise. 

3.1.  Field  amplitude  and  phase  on  the  pupil  plane 

We  start  by  establishing  a  relationship  between  the  angular  distribution  of  the  field  derived  in 
Eq.  (1)  and  a  geometric  ray  interpretation  to  demonstrate  the  distinct  characteristics  of 
subcritical  and  supercritical  angular  components  of  the  collected  light.  Each  angular  spectrum 
component  is  associated  with  a  light  ray  traveling  perpendicular  to  the  phase  front  of  the 
plane-wave  in  the  silicon  immersion  medium  [16,  29].  Figure  2  shows  the  intensity  profile 
across  the  pupil  plane  (Dp  =  2fohj NAobJ  =  5.56mm  where  Dp  refers  to  the  pupil  diameter)  for 

axial  dipole  locations  of  d  =  0,  d  =  2inJ2  and  d  =  llns  shown  left  to  right.  The  intensity  profile 
depends  strongly  on  the  axial  location  for  both  vertical  and  horizontal  dipoles  in  the  forbidden 
light  zone  (outside  the  black  circle)  unlike  the  allowed  light  zone  (inside  the  black  circle) 
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[30].  The  longitudinal  component  of  the  wave-vector  in  the  insulating  medium  ( k„im )  acquires 
imaginary  values  at  supercritical  angles  ( daSIL  >  sin  ( nms  / nSi ) )  and  therefore  the  amplitude 
of  the  supercritical  plane  wave  components  in  the  immersion  medium  decays  rapidly  as  a 
function  of  polar  angle  6obJ  (see  Figs.  2(d)  and  2(e)). 


Fig.  2.  The  logarithm  of  the  electric  field  intensity  map  (ignoring  the  constant  crpo)  is  shown 
for  horizontal  dipole  (left  column)  and  vertical  dipole  (middle  column)  when  d  =  0  (a),  d  = 

XinJ2  (b),  d  =  Xins.  The  black  rings  inside  the  images  denote  the  circle  corresponding  to  the 
critical  angle.  The  region  inside  and  outside  of  the  ring  corresponds  to  the  allowed  and 
forbidden  light  zones,  respectively.  The  size  of  each  image  is  8  mm  by  8  mm.  (d)  A  cross 

section  of  the  intensity  profile  as  a  function  of  polar  angle  (  Qoh.  )  at  q>  =  0  for  horizontal  dipole. 

(e)  The  same  as  in  (d)  except  for  vertical  dipole. 

The  evanescent  wave  origin  of  the  propagating  waves  at  the  supercritical  angles  impacts 
the  phase  front  of  the  collected  light.  The  phase  of  a  pencil  of  light  at  a  given  polar  angle  at 

the  pupil  is  defined:  'Ppupil(0)  =  tan-1  (  Im(Ecffx)/ Re(EGK5)j .  The  wave-front  aberration 
originated  from  the  forbidden  light  is  found  by  subtracting  the  spherical  aberration  term 
from  the  total  phase:  'FIL  =  "V pllpli  - 'F,ph  [31,  32].  For  an  aberration-free 

optical  system,  *Fpupil  is  expected  to  be  a  constant  and  thus  independent  of  the  polar  angle. 

Figure  3  plots  'Frl  as  a  function  of  polar  angle  greater  than  the  critical  angle  of  TIR  on  the 
Gaussian  reference  sphere  ( Oob. )  and  dipole  distance  from  the  dielectric  interface  (d).  Even  if 

the  imaging  system  is  free  of  optical  aberrations,  collection  of  evanescent  waves  introduces 
~ltnJ5  peak-to-valley  wave-front  distortion  at  high  aperture  angles,  degrading  the  system 
performance.  Such  phase  distortion  at  supercritical  angles  originates  from  the  complex  valued 
Fresnel  transmission  coefficients,  a  condition  similarly  observed  in  Goos-Hanchen  effect  that 

explains  the  relative  phase  shift  of  a  beam  undergoing  TIR  [33].  Note  that  'Fh  is  independent 
of  the  dipole  height  as  the  phase  distortion  originates  from  the  complex  Fresnel  coefficients 
( t'"s  and  f“  )  in  O  terms  in  Eq.  (3).  Furthermore,  the  phase  discontinuity  in  term  at  the 

critical  angle  leads  EGRS  component  of  the  horizontal  dipole  to  have  an  abrupt  n/2  phase  shift. 
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(a)  (b)  (c) 


Fig.  3.  Wave-front  aberration  introduced  by  the  forbidden  light  (  t  fl  in  units  of  radian)  in  (a) 

EgRS  ,  (b)  E™S  of  horizontal  dipole  and  (c)  E^RS  of  the  vertical  dipole  as  a  function  of 
polar  angle  0obj  and  the  distance  d  in  units  of  waves.  Note  that  vertical  dipole  does  not  have 
the  E°RS  component. 

3.2.  Image  of  dipoles  on  the  detector 

We  study  the  optical  images  of  dipolar  objects  on  a  wide-field  detector  using  the  formalism 
developed  in  Theory  section.  The  first  row  in  Fig.  4  shows  images  of  a  horizontal  dipole 
located  on  the  interface  when  all  light  in  subcritical  and  supercritical  angles  is  collected  (a), 
only  allowed  light  is  collected  (b)  and  only  forbidden  light  is  collected  (c).  The  peak  intensity 
on  the  detector  originating  from  the  forbidden  light  is  approximately  an  order  of  magnitude 
greater  as  compared  to  the  allowed  light.  The  overall  response  in  (a)  is  thus  dominated  by  the 
forbidden  light  leading  to  a  spot  size  of  approximately  kinJA.  For  a  horizontal  dipole  located  at 
a  distance  of  d  =  ).ins  from  the  interface,  the  peak  intensity  of  the  image  originating  from  the 
forbidden  light  is  approximately  two  orders  of  magnitude  smaller  as  compared  to  the  allowed 
light  as  shown  in  the  second  row  in  Fig.  4.  The  drop  in  the  collected  forbidden  light  in  the 
supercritical  aperture  angles  (compare  Figs.  2(a)  and  2(c))  as  well  as  the  spherical  aberration 
leads  to  reduction  in  image  brightness  and  spot  confinement  even  if  the  dipole  is  in  the 
geometric  focus  at  the  aplanatic  point. 


(d)  (e)  (0 


Fig.  4.  Normalized  wide-field  detector  images  of  a  horizontal  dipole  when  d  =  0  (top  row)  and 
d  =  A inS  (bottom  row).  Intensity  images  for:  (a)  and  (d)  full  collection  NA  of  the  microscope,  (b) 
and  (e)  only  subcritical  angle  components  (allowed  light),  (c)  and  (f)  only  supercritical  angle 
components  (forbidden  light).  The  edge  length  of  each  image  is  3 Ains  x  magnification.  The 
color  scales  are  normalized  to  the  maximum  intensity  in  (a). 
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Fig.  5.  Same  as  in  Fig.  4  except  for  a  vertical  dipole. 

We  repeat  the  simulations  for  a  vertical  dipole  and  the  results  are  shown  in  Fig.  5.  The 
vertical  dipole  images  similarly  follow  the  conclusions  obtained  for  horizontal  dipole 
regarding  the  effect  of  the  forbidden  light.  When  the  vertical  dipole  is  located  at  the  interface, 
the  contribution  from  the  forbidden  light  dominates  the  overall  image  and  thus  leads  to  a 
confined  spot  geometry.  At  a  distance  of  d  =  ?.ins,  the  spot  enlarges  significantly  as  the 
contribution  of  forbidden  light  becomes  negligible  compared  to  the  allowed  light 

Figure  6  plots  the  spot  size  and  relative  collection  efficiency  as  a  function  of  dipole  height. 
The  intensity  image  of  the  dipole  is  integrated  on  the  detector  to  calculate  the  metric  of 
relative  collection  efficiency  and  the  spot  size  is  defined  as  the  full-width-at-half-maximum  of 
the  peak  and  the  dip  for  the  horizontal  and  vertical  dipoles,  respectively.  The  collection  of  the 
forbidden  light  at  supercritical  aperture  angles  leads  to  a  spot  size  as  small  as  0.21,„s  and 
0. 1 8/',„s  for  horizontal  and  vertical  dipoles,  respectively.  The  contribution  of  the  forbidden 
light  becomes  negligible  with  respect  to  the  allowed  light  when  the  dipole  height  reaches 
~0.5kjns  leading  to  approximately  a  factor  of  5  drop  in  collection  efficiency  as  well  as  loss  of 
the  sub-diffraction  spatial  resolution.  Note  that  the  circular  symmetry  of  the  spot  depends  on 
the  distance  of  the  horizontal  dipole  from  the  interface  unlike  the  vertical  dipole.  The  circular 
symmetry  of  the  spot  on  the  detector  is  preserved  for  dipole  height  smaller  than  ~0.5/'„„.,  it 
however  transforms  into  an  elliptical  spot  for  larger  values  of  dipole  height. 
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Fig.  6.  Collection  efficiency  shown  in  intensity  axis  and  spot  size  (in  x  and  y  axes)  as  a 
function  of  dipole  height  d  for  horizontal  (left)  and  vertical  (right)  dipoles.  In  the  vertical 
dipole  case,  the  spot  geometry  is  circularly  symmetric  therefore  the  lines  showing  the  spot  size 
in  x  and  y  axes  overlap.  X  refers  to  the  wavelength  in  insulating  material.  FL  and  AL  refer  to 
forbidden  light  and  allowed  light. 

Two-point  resolution  refers  to  the  ability  of  an  optical  system  to  resolve  two  adjacent 
objects  [34].  Such  optical  analysis  is  practically  important  for  fault  localization  through 
distinguishing  the  defective  feature  from  a  neighboring  one  in  IC  failure  analysis  applications. 
Figure  7  shows  the  images  of  two  non-interacting  horizontal  dipoles  separated  by  a  distance 
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of  0.351,„j  in  the  limiting  cases  of  complete  coherent  (in-phase)  and  incoherent  dipoles. 
According  to  the  Sparrow  Criterion,  the  dipolar  objects  are  resolved  both  in  coherent  and 
incoherent  cases  when  they  are  both  located  on  the  dielectric  interface  ( d  =  0)  as  shown  in 
Fig.  7(a)  [35].  A  particularly  interesting  case  is  when  one  of  the  dipole  is  located  at  a  different 
depth  than  the  other  dipole.  Figure  7(b)  illustrates  this  condition  when  the  dipole  on  the  left  is 
located  at  d  =  87.5nm  (0. 12, „v)  (typical  separation  between  adjacent  interconnect  levels)  below 
the  right  dipole  situated  on  the  interface.  The  reduction  in  the  forbidden  light  collected  from 
the  left  dipole  renders  the  peaks  no  longer  resolved  in  the  coherent  case.  Despite  the  loss  of 
the  modulation  depth,  the  resolution  is  retained  in  the  incoherent  condition  and  the  dip 
between  the  two  peaks  shifted  towards  the  dipole  on  the  left  [36].  When  both  dipoles  are 
situated  at  d=  87nm  (0. 1 they  become  completely  unresolved  in  the  coherent  case  and  the 
modulation  depth  between  the  peaks  of  the  incoherent  dipoles  is  reduced  by  -30%. 


Fig.  7.  Normalized  wide-field  detector  images  of  two  dipoles  separated  by  0.35/1;, „  when  (a) 
both  dipoles  positioned  at  d  =  0,  (b)  dipole  on  the  right  at  d  =  0  and  the  dipole  on  the  left  at 
d  =  0.12IWS.  First  column  corresponds  to  coherent  case,  middle  column  corresponds  to 
incoherent  case  and  the  third  column  is  for  the  cross  sections  of  the  two  conditions  along  the 
lateral  axis.  The  edge  length  of  each  image  is  2Ains  x  magnification. 


3.3.  Image  of  two  dimensional  objects  on  the  detector 

We  simulate  incoherent  wide-field  images  of  two-dimensional  objects  buried  in  the  insulating 
medium  to  study  the  imaging  performance  on  test  objects  encountered  in  IC  chips.  We 
assume  the  horizontally  polarized  incoherent  illumination  is  uniform  in  the  field  of  view  and 
the  optical  system  has  the  configuration  explained  at  the  beginning  of  the  section. 

A  typical  test  object  for  the  optical  analysis  is  critical  dimension  features  for  estimating 
the  spatial  resolution  and  contrast  at  various  interconnect  levels.  Figure  8  shows  the  optical 
image  of  such  a  test  structure  located  at  an  axial  distance  of  d  =  0,  438nm  (0.5/.,„v)  and  876nm 
(Ims)  from  the  dielectric  interface,  from  left  to  right,  respectively.  The  line  pitch  ranges  from 
175nm  (0.22,„s)  to  700nm  (0.81,„s)  with  an  increment  of  87.5nm  (0.1  ).ins)  between  each  group. 
All  groups  except  the  first  one  are  well  resolved  according  to  Sparrow  criterion  i n  d  =  0  case. 
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The  smallest  groups  resolved  in  the  latter  cases  are  the  fourth  and  the  seventh  groups, 
respectively.  Note  that  the  difference  of  image  brightness  between  Figs.  8(a)  and  8(b)  is  more 
drastic  as  compared  to  the  difference  between  Figs.  8(b)  and  8(c)  due  to  the  exponential  decay 
of  the  forbidden  light  intensity  as  a  function  of  the  object  distance  from  the  interface. 


Fig.  8.  Simulated  images  of  a  resolution  target  with  periodic  grating  structures  at  different 
depth  with  respect  to  the  dielectric  interface,  (a)  d  =  0,  (b)  d  =  438nm  (XinJ2),  (c)  d  =  876nm 
The  size  of  each  image  is  7.32pm  x  15.52pm  (8.3/L,n5  x  17.72/WJ)  x  magnification. 

In  practice,  closely  spaced  objects  are  typically  found  at  different  interconnect  levels  in 
ICs.  Figure  9(a)  illustrates  an  example  of  metal  wiring  with  a  constant  line  pitch  of  876nm 
(Xins)  in  the  transverse  plane  and  located  at  d  =  438nm  ().iml2)  and  d  =  876nm  (2,„s),  shown  in 
red  and  blue  colors,  respectively.  Figure  9(b)  shows  the  simulated  optical  image  of  the  metal 
wiring  assuming  the  aplanatic  plane  coincides  with  the  plane  of  LI  and  L2.  According  to  the 
vertical  linecut  obtained  from  the  optical  image,  LI  and  L2  are  well  resolved  according  to 
Sparrow’s  criterion  at  which  the  valley  between  the  peaks  reaching  down  to  -0.45  whereas 
the  modulation  contrast  between  the  peaks  of  L3  and  L4  is  almost  vanishing  (-0.02).  The 
slight  difference  in  the  intensity  level  of  the  dip  between  L2-L3  and  L3-L4  peaks  originates 
from  the  unequal  contrast  level  of  the  L2  wire  with  respect  to  L3  and  L4,  as  explained  in  the 
two-point  resolution  case  study. 


LI  L2 


(c) 


Fig.  9.  Layout  and  the  simulated  images  of  a  buried  two-level  metal  wiring,  (a)  The  red  and 
blue  lines  are  located  at  a  distance  of  XirJ2  and  Xins,  respectively.  The  scalebar  corresponds  to  a 
length  of  876  nm  (AJns).  (b)  The  optical  image  of  the  object,  (c)  Cross  section  of  the  image 
along  the  dotted  line  shown  in  (b).  The  size  of  the  images  corresponds  to  a  field  of  view  of 
8.86pm  x  6.38pm  (10.1  \A.ins  x  7.282/W5)  x  magnification. 
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The  examples  above  are  chosen  considering  the  dielectric  structure  of  the  conventional 
bulk  silicon  process  technology  in  which  the  interconnect  layers  are  buried  beyond  a  single 
dielectric  interface  between  the  silicon  substrate  and  the  insulating  media.  Our  model  can  be 
also  applied  to  more  complex  layered  media  applications  such  as  silicon-on-insulator  chips  in 
which  the  silicon  junction  is  contained  between  a  thin  buried  insulating  box  and  the  insulating 
media  surrounding  the  interconnect  layers.  In  order  to  account  for  the  multiple  reflections 
occurring  in  the  stratified  media,  the  generalized  Fresnel  coefficients  can  be  used  in  place  of 
the  Fresnel  coefficients  used  here  [37].  Besides  the  concerns  originating  from  the  evanescent 
modes,  practical  complications  further  affect  system  performance  such  as  the  optical 
aberrations  originating  from  imperfect  solid  immersion  lens  dimensions,  spectral  bandwidth 
of  the  light  and  the  change  in  the  operating  temperature  of  the  chips.  Such  practical  issues  will 
be  investigated  in  future  through  adapting  the  Green’s  function  model  developed  in  this  study. 

4.  Conclusion 

In  this  theoretical  study  of  forbidden  light  in  high  NA  aSIL  microscopy,  we  demonstrated  that 
the  dielectric  interface  between  the  silicon  substrate  and  insulating  medium  surrounding  the 
circuit  features  has  a  major  impact  on  the  imaging  performance.  Circuit  features  located 
within  a  half-wavelength  of  the  interface,  comprising  the  gate  and  first  few  metal  levels  of  the 
modern  IC  chips,  can  be  interrogated  effectively  with  sub-diffraction  transverse  spatial 
resolution  and  high  photon  collection  efficiency  due  to  the  forbidden  light  collected  by  the 
aSIL  microscope.  For  objects  buried  farther  into  the  insulating  medium,  only  light  from  the 
allowed  light  zone  can  be  collected,  leading  to  up  to  a  factor  of  5  drop  in  light  collection 
efficiency.  In  addition,  the  diffraction  limited  spatial  resolution  is  compromised  due  to  optical 
aberrations  imposed  by  the  dielectric  interface. 

The  developed  vector  field  theory  for  evaluating  the  effects  of  evanescent  waves  in  aSIL 
microscopy  of  ICs  can  assist  in  quantitative  analysis  of  optical  images  in  the  context  of  fault 
isolation  and  failure  analysis  as  well  as  designing  adaptive  optics  for  compensating  optical 
aberrations  for  high  resolution  optical  inspection  as  pure  ray-tracing  model  based  software 
cannot  account  for  the  behavior  of  the  evanescent  waves.  Even  though  the  focus  of  this  study 
is  limited  to  semiconductor  IC  applications  in  which  the  resonant  scattering  is  usually  of 
interest,  the  model  can  be  adapted  to  applications  in  fields  such  as  quantum  optics  and 
biophotonics  for  which  the  non-resonant  scattering  processes  are  also  important. 
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